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ABSTRACT 

We report on the first X-ray observations of the neutron star soft X-ray transient (SXT) XTE J2123-058 in 
quiescence, made by the Chandra X-ray Observatory and BeppoSAX, as well as contemporaneous optical 
observations. In 2002, the Chandra spectrum of XTE J2123-058 is consistent with a power-law model, or 
the combination of a blackbody plus a power-law, but it is not well-described by a pure blackbody. Using the 
interstellar value of A^h, the power-law fit gives F = 3.1^o^ and indicates a 0.3-8 keV unabsorbed luminosity 
of (9;!^) X 10^' (ii/8.5 kpc)^ ergs s"' (90% confidence errors). Fits with models consisting of thermal plus 
power-law components indicate that the upper limit on the temperature of a 1.4M0, 10 km radius neutron star 
with a hydrogen atmosphere is kT^n < 66 eV, and the upper limit on the unabsorbed, bolometric luminosity is 
Loo < 1 .4 X lO"*^ ergs s"' , assuming li = 8.5 kpc. Of the neutron star SXTs that exhibit short (< 1 year) outbursts, 
including Aql X-1, 4U 1608-522, Cen X-4, and SAX J18 10.8-2609, the lowest temperatures and luminosities 
are found for XTE J2123-058 and SAX J1810.8-2609. From the BeppoSAX observation of XTE J2123-058 
in 2000, we obtained an upper limit on the 1-10 keV unabsorbed luminosity of 9 x 10^^ ergs s~'. Although 
this upper limit allows that the X-ray luminosity may have decreased between 2000 and 2002, that possibility 
is not supported by our contemporaneous /?-band observations, which indicate that the optical flux increased 
significantly. Motivated by the theory of deep crustal heating by Brown and co-workers, we characterize the 
outburst histories of the five SXTs. The low quiescent luminosity for XTE J2123-058 is consistent with the 
theory of deep crustal heating without requiring enhanced neutron star cooling if the outburst recurrence time 
is ^ 70 years. 

Subject headings: accretion, accretion disks — stars: individual (XTE J2123-058, SAX J1810. 8-2609) — 
stars: neutron — X-rays: stars 



1. INTRODUCTION 

Accreting neutron stars can be found in high-mass (HMXB) 
or low-mass (LMXB) X-ray binary systems. The majority of 
HMXBs have transient X-ray emission. Their outburst spec- 
tra are relatively hard and X-ray pulsations from these highly 
magnetized (B 10'^ G) neutron stars are typically detected. 
A wide variety of X-ray behaviors are seen for neutron star 
LMXBs, but, in general, the lack of X-ray pulsations from 
most (but not all) of these systems, and the emission of type I 
X-ray bursts from some, suggest that they harbor neutron stars 
with relatively low magnetic field strengths {B ^ 10**"^ G). 
During w 33 years of X-ray observations, some sources (e.g., 
Sco X-1, Cyg X-2) have maintained luminosities approaching 
the Eddington limit of « 10^'^ ergs s"', while others are able 
to maintain persistent luminosities several orders of magni- 
tude lower (Wilson et al. 2003). In addition, there is a class of 
transient neutron star LMXBs for which the luminosity varies 
from a substantial fraction of Eddington to quiescent levels 
typically near lO-'^"-''' ergs s"' . In outburst, these systems have 
relatively soft spectra compared to the HMXBs, and are com- 
monly grouped with black hole transients as soft X-ray tran- 
sients (SXTs). 

In quiescence, most neutron star SXTs exhibit X-ray en- 
ergy spectra with a component that is typically fitted well 
by a blackbody, suggesting that the origin of this compo- 
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nent is thermal emission from the surface of a cooling neu- 
tron star Although a pure blackbody often provides a good 
fit to the spectrum, unphysical neutron star radii near 1 km 
are in ferred unless an atmosphere is modeled ( Rutledge et aO 
Il999h . In addition to the thermal component, the energy 
spectra often contain a second component that has a power- 
law shape. The brightest and best studied systems in this 
class, Cen X-4 and Aql X-1, usually display both compo- 
nentejRutledae et al. 2001, 2002a; Campana & Stefla 2003J 
ICampana et al. 2004). However, other systems may be dom- 
inated by the thermal component, such as MXB 1659-29 
( IWiinands e t al. 2003b) and sources X-5 and X-7 in the glob- 
ular cluster 47 Tucanae (Heinke et al. 2003), or by the power- 
law c omponent, such as SAX J1808. 4-3658 ( Ca mpana et alJ 
I2OOI and EXO 1745-248 ( Wiinands et al. 2003a). Although 
theories for the t hermal co mponent, such as the deep crustal 
heating model ofBro wnTBildsten & Rutledga(il998i) . are rel- 
atively well-developed and are being tested with observa- 
tions, the origin of the power-law component is not under- 
stood beyond suggestions that it may be related to accretion 
onto the neutron star magnetosphere (Campana et al. 199^ 
or a putative pulsar wind colliding with infalling matter from 
the companion star (Tavani 1991). In addition to our lack 
of understanding of the power-law component, questions re- 
main about the mass accretion rate in quiescence, the ori- 
gin of rapid (100-10, 000 s) variabiUty (Ikutledge et al. 2002^ 
ICampana et alJ2004l) . and the origin of variability in the ther- 
mal component on longer time scales (Rutledge et al. 2002a). 
Another important question is if quiescent observational prop- 
erties correlate with other known differences between neutron 
star SXTs, such as whether they are millisecond X-ray pulsars 
(during outbursts) or not, whether the systems are in the field 
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or in globular clusters, and whether their X-ray outbursts are 
long (years to decades) or short (weeks to months). 

Here, we report on X-ray and optical observations of the 
field neutron star SXT XTE J2123-058 taken during qui- 
escence. XTE J2123-058 ha d its only detected X-ray out- 
burst in 1998 June- August (iLevine. Swank & S mith 19981 
iTomsick etal] ri999). and we focus on observations made 
with Chandra, BeppoSAX, and optical telescopes 2-A years 
after the outburst. During the outburst, the Rossi X-ray 
Timing Explorer (RXTE) detected type I X-ray bursts and a 
pair of kHz quasi-periodic oscillations (Homanetal. 1999; 
l^msicketal. 1999), indicating that the system contains 
a rapidly rotating neutron star However, coherent X-ray 
pulsations were not found. The 6 hr binary orbital pe- 
riod and the fact that the binary inclination of the sys- 
tem is relatively high were established from optical mod- 
ulation and the presence of partia l echpses in the opti- 
cal light curve (Tomsick et al. 'ri9 99HSOTir Wu & Gallowavl 
[1999; Zuritaetal. 2000; Shahbaz etall I2003h . The high 
Galactic latitude {b = -36°) and low extinction have allowed 
for detailed optical studies of XTE J2123-058 in quiescence 
even though the source is rather faint at its relatively large 
distance (8.5 ±2.5 kpc). The optical observations show that 
XTE J2123-058 consists of a K7 V star on or close to the 
main sequence and a neutron star for which mass determi- 
nations of 1.5 ±0.3 M(7) (Tomsick et al. 2001; Casares et aP 
12002; Tomsick et al. 2002) and 1.04-1. 56 ( Shahb az etaD 
12003.) have been obtained. The focus of this paper is the first 
X-ray study of XTE J2123-058 in quiescence. 

2. OBSERVATIONS AND ANALYSIS 

We observed XTE J2123-058 with Chandra on UT 2002 
November 13 (ObsID 2709), using the Advanced CCD Imag- 
ing Spectrometer (ACIS) with the target placed on one of 
the back-illuminated ACIS chips (ACIS-S3). For our anal- 
ysis, we used the "level 2" event list produced by the stan- 
dard data processing with ASCDS version 6.9.2 using Cal- 
ibration Data Base (CALDB) version 2.17. Light curves 
using counts from the full field-of-view do not show any 
background flares, allowing us to use the data from the full 
17,706 s integration. Using the Chandra Interactive Analysis 
of Ob servations (CIAO) v ersion 3.0 software routine wavde- 
tect tPreeman et al.l2002l) . we searched for sources on the S3 
chip in the 0.3-8 keV energy band. We detected 22 sources 
with counts between 5 and 108 per source, using a detec- 
tion threshold of lO"*", including a 24 count source at R.A. 
= 21*'23"14154, decl. = -05°47'53".2 (equinox J2000, un- 
certainty 0".6). This position i s con sistent with the target's 
optical position ( Tomsick et al.1 ll999h . and we conclude that 
this source is the quiescent X-ray counterpart of XTE J2123- 
058. 

We also report on a BeppoSAX observation of the XTE 
J2123-058 field made on UT 2000 May 12-13. We produced 
1-10 keV images using the data from the two Medium En- 
ergy Concentrator/Spectrometer (MECS) units that were op- 
erational during the observation (units 2 and 3). We also pro- 
duced a 0.1-10 keV image using data from the Low Energy 
Concentrator/Spectrometer (LECS). We obtained a MECS 
exposure time of 46,340 s and a LECS exposure time of 
17,080 s. To search for sources, we convolved each of the 
three images with a two-dimensional Gaussian with a width 
(cr) of 5 pixels (40") in both directions. Only one source 
was clearly detected in both of the MECS images at R.A. = 
2l''22"50^6, decl. =-05°45'09" (equinox J2000, uncertainty 



- 1'), and this source, which we call SAX J2122. 8-0575, is 
also present in the LECS image. It is clear that SAX J2122.8- 
0575 is not XTE J2123-058 as they are separated by 6'.6. We 
conclude that XTE J2123-058 was not detected during the 
BeppoSAX observation, and we derive an upper limit on its 
X-ray flux below. 

Finally, we obtained R-hmd images on three occasions 
close to the times of the X-ray observations. As shown in 
TablelU we observed XTE J2123-058 using the 2.4 m Hiltner 
telescope of the MDM Observatory on 2000 July 24, about 
2 months after the BeppoSAX observation. We also observed 
XTE J2123-058 with the Shane 3 m telescope of Lick Ob- 
servatory about 2 months before the Chandra observation and 
again at MDM about 2 weeks after the Chandra observation. 
For both MDM runs, we used the same SITe 2048 x 2048 
pixel, thinned, back-illuminated CCD with a spatial scale of 
0".275 per 24 /im pixel, and an R filter that is very close to 
Harris R. '* At Lick, we used the Prime Focus Camera, with 
a SITe 2048 x 2048 pixel, thinned CCD with a spatial scale 
of 0".296 per 24 /im pixel, and a Kron-Cousins R filter. In 
all, we obtained thirteen 600-700 s exposures (see Table 0, 
and we reduced the images using standard IRAF ^ routines. 
For XTE J2123-058, we carried out the photometry with the 
IRAF package phot and used two calibrated reference stars 
with R magnitudes of 19.47 and 19.51. In 2002, the con- 
ditions at Lick were photometric, and we observed Lando^ 
(ir992) standards to obtain the calibration. We note that this 
calibration is about 0.1 magnitudes brighter than the calibra- 
tion previo usly obtained using M DM observations from 1998 
reported in lTomsick et alJ lll999ii) . 

3. ENERGY SPECTRUM AND SOURCE LUMINOSITY 

We used the CIAO software routine psextract to produce 
the ACIS energy spectrum for XTE J2123-058 and to create 
the appropriate instrument response matrix for the spectrum. 
The software used CALDB 2.23 to create the response matrix, 
and we included a correction for the time-dependence of the 
ACIS response. The source spectrum included counts from a 
circular region with a 5 pixel (about 2".5) radius, and we es- 
timated the background level using counts from a source-free 
annulus around the target position. The source spectrum con- 
sists of 24 counts, and we estimate a background level of 0.7 
counts in the extraction region. We produced a "light curve" 
of the source in six time bins of « 3000 s. Each bin contains 
between 2 and 6 counts, which is consistent with a constant 
flux; however, the low count rate does not allow us to place 
tight constraints on the possible amplitude of variability. The 
small number of counts also indicates that statistics are 
not appropriate for spectral analysis, as the assumption of a 
Gaussian probability distribution in each spectral bin is not 
met. Thus, we carried out ou r spectral analysis by minimiz- 
ing the Cash statistic llCa^ ri979). which is appropriate in 
cases where the assumption of a Poisson probability distribu- 
tion in each spectral bin is valid. We fitted the spectrum using 
XSPEC 11.2. 

The spectra of other neutron star SXTs are typically well- 
described by a blackbody, a power-law, or both components 
with interstellar absorption. We began by fitting the spec- 

See http://www.astro. lsa.umich.edu/obs/mdni/technical/filtersl for the ex- 
act transmission curve. 

^ IRAF (Image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Fig. \.— ChandralAClS energy spectrum of XTE J2123-058. Top: folded 
through the detector response. Bottom: unfolded. In each panel, the solid 
line is a power-law fit to the data (using Cash statistics as described in the 
text) with Wh fixed at the interstellar value. The power-law photon index is 
3.1:*}J g, and the unabsorbed 0.3-8 keV luminosity is (9;*) X 10^^' ergs s"'. A 
blackbody model (dashed line), by itself, does not provide a good description 
of the specfium. 



trum for XTE J2123-058 with an absorbed power-law, and 
the results are given in Table|2] Once we found the best-fitted 
parameters by minimizing the Cash statistic, we determined 
the quality of the fit and the 90% confidence errors on the 
parameters by producing and fitting 10,000 simulated spec- 
tra. We used the best-fitted parameters from the fits to the 
actual data as input to the simulations. Using simulations to 
determine the parameter errors is necessary only because of 
the small number of counts in the spectrum. We tested our 
method by producing spectra with 2-3 times as many counts, 
and we found that the errors produced by calculating changes 
in the Cash statistic (the standard technique) match the val- 
ues we obtain using simulations. In addition to the power-law 
model, we used simulations to determine the errors on param- 
eter for all other fits presented in this work. For the power-law 
model, we obtain A^h = i'^ ■Q'^'^q) x 10^" cm"- for the column 
density and F = 3.1;!;Qg for the photon index. The fit using 
a blackbody model is significantly worse than the power-law 
model as indicated by the fact that for 97% of the simulated 
spectra we obtained a better Cash statistic than the one ob- 
tained when fitting the actual spectrum (compared to 69% for 
the power-law model). 

We refitted the spectrum with the same spectral mod- 
els but with A^H fixed to the interstellar column density. 
The interstellar A ^h comes from the Ay measurement of 
iHynes et al.' ("20011), which gives A^h = (6.6 ±2.7) x 10^° cm'^, 
and the total Galac t ic H I value of 5.7 x 10^° cm"^ from 
iDickev & Lockmaiil J1990I) . We adopt a value of 6 x 10^" 
cm"^. For the power- law alone, we obtain F = 3.1;!;Qg and an 
unabsorbed 0.3-8 keV flux of (l.l!|j;^) x lO"'"* ergs cm'^ s"'. 



which corresponds to a luminosity of (9;I;3) x 10^' ergs s"' 
at a distance of 8.5 kpc. The best fitted blackbody temper- 
ature is 240;!;5Q '■^^ quality of the blackbody fit is even 
worse (98.7%) with A^h fixed. Figure^shows the spectrum 
(rebinned for clarity, although the Cash fits are performed 
without rebinning) along with the best fitted power-law and 
blackbody models with A^h fixed. It is clear that the curvature 
of the blackbody model is too large, underpredicting the ob- 
served spectrum below 0.7 keV and above 2 keV. This, along 
with the measurements of the fit quality from the simulations, 
indicates that the spectrum is well-described by a power-law 
model or the combination of a blackbody and a power-law 
model, but it is not well-described by a blackbody alone. 

Although a power-law alone provides an adequate descrip- 
tion of the ACIS spectrum, we performed additional fits to 
obtain limits on the temperature and luminosity of the puta- 
tive thermal component that is expected to be emitted from 
the surface of the neutron star. We fitted the spectrum with a 
model consisting of a power-law component along with ther- 
mal emission from an atmosphere composed of hydrogen. 
The latter component was model ed usin g the "Neutron Star 
Atmosphere" (NSA) model of IZavlin. Pavlov & ShibanoJ 
in"996t) . We assumed a neutron star mass of 1.4M0, which 
is consistent with the measured mass, and a radius of 10 km. 
With these parameters fixed, the remaining free parameters in 
the NSA model are the temperature {kT^fi) and the distance. 
For XTE J2123-058, there are several a rguments that lead to 
distan ce estimate s in the range 5 -15 kp c JTomsick et al.ll999t 
Homa n et al.l[T99 9: Zurita et al] l200fll) . but flie most rehable 
estimates come from optical observations of the source in qui- 
escence. Previous estimates include 8 ±3 kpc (Zurita et a^ 
120001). 8.5±2 5kpc dTomsick et al. 2001), and 9.6 ±1.3 kpc 
iCasares et al.l20(^ . Here, we adopt a range of 6-1 1 kpc. We 
performed the NSA plus power-law fits with the distance fixed 
to three values spanning this range (6, 8.5, and 11 kpc). As 
shown in Table|3] the 90% confidence upper limits on kT^ff are 
57, 66, and 73 eV, respectively, for these three distances. Us- 
ing Loo = '^TiR-aT^iid-lGM/Rc^), where « = 10 km and M 
= 1 .4 Mq, upper limits on the unabsorbed luminosity from the 
NSA component as seen by a distant observer are 8.0 x 10^', 
1.4 x lO-'^ and 2.1 x lO-'^ ergs s~' for distances of 6, 8.5, and 
1 1 kpc, respectively. While these values are bolometric lumi- 
nosities, it should be noted that we cannot rule out the pos- 
sibility that there is thermal emission at energies below the 
Chandra bandpass. 

Although XTE J2 123-05 8 was not detected by BeppoSAX 
in 2000 May, we can calculate upper limits on its luminos- 
ity if we assume that the energy spectrum was similar to that 
measured by Chandra. To calculate the upper limits, we as- 
sume the simplest model of an absorbed power-law with a 
photon index of 3.1 and the interstellar column density. For 
the MFCS (units 2 and 3 combined), we measure a 1-10 keV 
count rate of 1 .684 ks"' in a circle of radius 2' centered on the 
XTE J2 123-058 position. This is not significantly higher than 
the expected background rate (using a blank-sky pointing) of 
1.675 ks"'. These values imply a 3-(t upper limit on the count 
rate from XTE J2123-058 of 0.676 ks"', which corresponds 
to an unabsorbed flux < 1.02 x 10"'-' ergs cm"^ s"', and a 1- 
10 keV luminosity < 9 x 10^^ ergs s"' for = 8.5 kpc. For the 
LECS, the 0.1-1 keV count rate in a 2' radius circle centered 
on XTE J2123-058 is 1.230 ks"', which is actually somewhat 
lower than the expected background rate of 1 .466 ks"' . Using 
the expected background rate, the 3-(t upper limit on the un- 
absorbed 0.1-1 keV flux is 1.75 x 10"'^ ergs cm"^ s"', corre- 
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spending to a luminosity < 1 .5 x 10^"^ ergs s ' for d=S.5 kpc. 

4. OPTICAL RESULTS 

Figure|2|shows the R magnitudes for XTE J2123-058 from 
three exposures taken in 2000 July at MDM, two exposures 
taken in 2002 September at Lick, and eight exposures taken 
in 2002 November at MDM. Previously, Shahbaz et al. ( 2003) 
reported on extensive 7? -band photometry of XTE J2123-058 
taken between 1999 June and 2000 August. The /?-band 
light curve was relatively stable at that time, showing or- 
bital ellipsoidal modulations with a peak-to-peak amplitude 
of about 0.25 magnitudes. The dashed lines in Figure |2] in- 
dicate the range of the modulation in 1999-2000, which was 
21.70 <R< 21.95. While our 2000 May measurements of 
21.63 ±0.07, 21.75 ±0.08, and 21.57 ±0.07 are consistent 
with the brighter end of this range, our later observations indi- 
cate that XTE J2123-058 brightened between 2000 and 2002. 
The range of between 21.41 ±0.07 and 21.01 ±0.03 over 
25% of the 6 hr binary orbital period during the 2002 Novem- 
ber run indicates that the source continues to vary, but with 
a peak level that is at least 0.6-0.7 magnitudes brighte r tha n 
the peak level of /? = 2 1 .7 measured by Shah baz et al.l ( l2003l) . 
As our observations do not cover the full binary orbit, it is 
not clear whether this change is due to the addition of a con- 
stant component to the light curve or if the shape of the mod- 
ulations has changed. Although it would be useful to know 
the orbital phases corresponding to the 2002 data, the binary 
e phemeris is too uncerta in to extrapolate to 2002 as discussed 
in lTomsick et alJ (I2002h . In the future, it may be worthwhile 
to obtain an X-ray observation while measuring the optical or 
IR light curve for a full 6 hr binary orbit. 

As we are reporting the first sensitive X-ray observations 
of XTE J2123-058 in quiescence, comparison can be made 
with a prior prediction for the quiescent X-ray luminosity of 
this source that was made by modeling quiescent /?-band light 
curves in 1999 and 2000. Shahbaz et al. (2003) predicted an 
X-ray luminosity of ^lO"'^ ergs s"', which is only slightly 
higher than our upper limit from the BeppoSAX observation 
made in 2000, but an order of magnitude higher than we ob- 
served with Chandra in 2002. While the X-ray observations 
would allow for the possibility that the quiescent X-ray lu- 
minosity decreased significantly between 2000 and 2002, our 
optical observations made in 2000 and 2002 indicate that the 
7?-band flux actually increased over this time. If the 2000 R- 
band light curve showed a significant contribution from X- 
ray heating, it is difficult to see how a drop in X-ray flux 
could lead to an increase in the optical, thus, it is possible that 
ISha hbaz e t al. (.2003.) over-estimated the contribution from X- 
ray heating. 

5. INTERPRETATION 

X-ray observations of neutron star SXTs in quiescence 
provide tests of theoretical models for the thermal com- 
po nent. According to the theory of deep crustal heating 
by iBrown. Bildsten & Rutledg ^ (Il998t . the temperature of 
the neutron star core is maintained by nuclear reactions in 
the deep crust that occur when the mass accretion rate is 
high during outburst. A s the t hermal time scale for the 
core is --10,000 years ( Colpi et al.l 1200 ll) . the level of qui- 
escent thermal emission is set by the average mass accre- 
tion rate over this time span according to Lq = 9 x 10^^ 
<M>_ii ergs s"' (see Equation 1 of Rutledge et al. 2002b), 
assuming 1.45 MeV of heat deposited in the crust per ac- 
creted nucleon (iHaensel & Zdunikl[l990h . Here <M>-u is 
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Fig. 2.— R-band measurements ofXTEJ2123-058. (o) 2000 July at MDM 
Observatoiy. (h) 2002 September at Lick Observatory, (c) 2002 Novem- 
ber at MDM Observatory. The dashed lines delimit the levels measured by 
IShahbaz etali 112003) during 1999-2000. An increase in the R-band flux be- 
tween 2000 and 2002 is apparent. 



the mass accretion rate averaged over the thermal time scale 
of the core in units of 10~" M© yr"\ and Lq is the quiescent 
bolometric luminosity. While the time-averaged mass accre- 
tion rate is only a predictor of Lq if the neutron star core has 
reached thermal equilibrium, the lifetimes of LMXB systems 
are much longer than the thermal time scale of the core, and it 
is expected that thermal equilibrium has been established for 
all or nearly all of the LMXBs. While the long thermal time 
scale for the core precludes appreciable changes in the core 
temperature during a single SXT outburst, the neutron star 
crust can be heated significantly for the systems with longer 
outbursts (years to decades) so that the quiescent emission is 
determined by the ev olution of the physical conditions in the 
crust (Rutl edge et alJl2002b.,) . However, for the systems with 
shorter outbursts (weeks to months), the properties of the qui- 
escent thermal emission are primarily set by the conditions in 
the neutron star core (Brown, Bildsten & Rutledge 1998). 

For a system like XTE J2123-058 that has undergone one 
w 40 day outburst, the deep crustal heating theory implies 
that quiescent thermal emission provides information about 
the neutron star core. Here, we compare the quiescent prop- 
erties of XTE J2123-058 to similar systems. The compar- 
ison group includes field neutron star SXTs with outbursts 
lasting less than 1 year for which quiescent X-ray observa- 
tions have been reported. We only consider LMXB systems 
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from which X-ray bursts have been detected, proving that the 
accreting object is a neutron star. We restrict the compari- 
son group to field systems because the outburst histories for 
most transients in globular clusters are uncertain due to source 
confusion in instruments with angular resolution worse than 
Chandra. Thus, our comparison group consists of the neu- 
tron star SXTs Aql X-1, 4U 1608-522, Cen X-4, and SAX 
J1810. 8-2609. We also compare the quiescent properties of 
XTE J2 123-058 to those of the millisecond X-ray pulsai- SAX 
J1808. 8-3658 as the quiescent X-ray luminosities of the two 
sources are comparable. 

Table |3 compares the parameters of the thermal compo- 
nent for the five systems, including kTf^ff, R, Loo, and the 
gravitational redshift parameter, g = ^ \-2GM / Rc^, assum- 
ing that all the neutron star masses are 1.4M0. For the 
spectral fits to XTE J2123-058, we set 7? = 10 km, and 
iJonker. '\^n ands & van der Klis (2004) make the same as- 
sumption for SAX J1810.8-2609. However, for the other 
sources, we use the best fitted radii, and th e parameters are 
taken from Rutledge et al. ( 1999, 2001, 2002a). The thermal 
component is detected for Cen X-4, Aql X-1, and 4U 1608- 
522, and a power-law component is also required for Cen X- 
4 and for some of the observations of Aql X-1. For XTE 
J2123-058 and SAX J1810.8-2609, the thermal component 
is not statistically required, and we report upper limits on the 
thermal component parameters for these two systems. Also, 
for all the sources we assume the best current distance mea- 
surement. The temperature upper limits for XTE J2123-058 
(66 eV) and SAX J18 10.8-2609 (72 eV) are about a factor 
of 2.5 lower than the 4U 1608-522 measurement of 170 eV, 
and the luminosity upper limits for XTE J2123-058 and SAX 
J1810. 8-2609 are about a factor of 50 lower than the maxi- 
mum quiescent luminosity measured for Aql X-1. The tem- 
perature and luminosity of Cen X-4 are intermediate between 
the hottest and most luminous sources (Aql X-1 and 4U 1608- 
522) and the coolest and least luminous sources (XTE J2123- 
058 and SAX J1810.8-2609). 

Here, we characterize the outburst histories of these 
five sources to determine if they are related to the 
quiescent thermal properties as is expected in the 
Bro wn. Bildsten & Rutledge (1998) theory. From a re- 
view of SXT light c urves covering the time period between 
1969 and 1996 ( Chen. Shrader & Liviolll997) and the RXTE 
All-Sky Monitor (ASM) covering 1996 to 2003, 4U 1608- 
522 had 16 outbursts with a mean peak X-ray luminosity, 
Ipeak, of 2.5 X 10" (<i/3.6 kpc)2 ergs s"'. Aql X-1 had 21 
outbursts with Zpeak = 3.6 X lO" {dl5 kpc)^ ergs s"' . Cen X-4 
had two outbursts (in 1969 and 1979) with Zpeak = 5.6 x lO" 
(d/1.2 kpc)^ ergs s"', and XTE J2123-058 had one outburst 
(in 1998) with Ipeak = 1.7 x lO" (d/9,.5 kpc)- ergs s"'. Like 
XTE J2123-058, SAX J1810.8-2609 had only one outburst, 
also in 1998. 

SAX J18 10.8-2609 was discovered by Nataluc ci etsU 
ll2000) : and the only reported X-ray detections of this source 
in outburst were fmm BeppoSAX and ROSA T in 1998 March 
dNatalucci etaP 120001 fOreiner et al."1999^. However, the 
ASM light curve shown in Figure |3j; indicates that SAX 
J1810.8-2609 was in outburst by 1998 January 4 (MJD 
50,817) and probably as early as 1997 December 18 (MJD 
50,800). The dashed lines in Figure |3j; mark the times 
of the BeppoSAX discovery at MJD 50,882-50,885 and the 
ROSAT observation at MJD 50,896, while it is clear that 
the source was considerably brighter prior to this. From 
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Fig. 3. — All-Sky Monitor (ASM) light curves for the four neutron star 
SXTs in our comparison group that had outbursts during the RXTE lifetime. 
Each point represents the average ASM count rate over 1 day. One Crab 
flux equals 75 ASM counts s"'. For SAX J1810.8-2609 (c), the vertical 
dashed lines mark the times of previously reported BeppoSAX and ROSAT 
observations. Aql X-1 and 4U 1608-522 have had many outbursts during the 
lifetime of RXTE, and the light curves in (a) and {b) are representative. 



the 1.5-12 keV ASM light curve, Lpeak = 1.0 x lO" (<i/4.9 
kpc)^ ergs s"^ Although it has been suggested that SAX 
J1810. 8-2609 had an unusually low outburst luminosity 
( iJonker. Wiinands & van der Klisil200l . the value we derive 
from the ASM measurements is comparable to the other neu- 
tron star SXTs. Figure |3 also shows the ASM light curve 
for the XTE J2123-058 outburst as well as sample outburst 
light curves for Aql X-1 and 4U 1608-522 (Cen X-4 has not 
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had an outburst during the KXTE Ufetime). These light curves 
demonstrate that individual outbursts from different sources 
are similar in duration (typically 40-80 days) and overall 
shape, altho ugh it should be noted that exceptional outbursts 
do occur (Chen. Shr ader & Livioll997l) . 

For each of the neutron star SXTs, the time-averaged mass 
accretion rate over the past 33 years can be expressed as 



<M>= 



^peak^^outburstX 

ec2(33yr) 



iLpeakA^ 



(1) 



where is the number of outbursts, foutbmst is the typical du- 
ration of an outburst from the source, / is a factor with a 
value less than 1 .0 that accounts for the shape of the outburst 
light curve, e is the fraction of the rest mass energy released 
from accreted matter, and s = foutbui-st//ec^(33yr). As the du- 
rations and light curve shapes are similar for the sources in 
our comparison group, s is approximately the same for these 
sources, and the time-averaged mass accretion rate over the 
past 33 years is proportional to LpeakA^- We estimate s us- 
ing Wburst = 60 days and e = 0.2, which is the value of e 
that is typically assumed for accretion onto a neutron star 
(|Rutledge et al. 2002b). Although the precise value of / de- 
pends on the exact shape of the outburst light curve, it can be 
approximated as the mean outburst flux divided by the peak 
flux. For the XTE J2123-058 ASM Hght curve, the mean flux 
over a 60 day period that includes the 1998 outburst is 41% 
of the peak level, and we use / = 0.4, giving .? = 1 . 1 x 10"^^ s^ 
cm~^. Using this value of s for all five sources and the values 
for Zpeak and the number of outbursts (A^) given above, we cal- 
culate estimates of the time-averaged mass accretion rate over 
the past 33 years. The values for <M> are given in Table|3 
and they show that this quantity is a relatively good predictor 
of Loo- Aql X-1 and 4U 1608-522 have the the highest values 
of both <M> and Loo, XTE J2123-058 and SAX J1810.8- 
2609 have the lowest, and Cen X-4 is intermediate. The ratio 
of the value of <M> for Aql X-1 to that for XTE J2123- 
058 is 45, which is consistent with their ratio of Lq^, which 
is > 38. Similarly, the Aql X-1 to SAX J1810.8-2609 ratio 
of <M> is 76, and this is consistent with their ratio of Loo, 
which is > 27. While these m easurements are in line with the 
iBrown. Bildsten & Rutled ge ( 1998) theory, a caveat is that we 
cannot be certain that X-ray flux history over the past 33 years 
reflects the behav ior over the last 10 ,000 years. 

We can use the lBrown. Bildsten '& Rutledg3 (ll998h theory 
to predict the recurrence time for outbursts of XTE J2123- 
058, assuming that the outbursts are similar. Using L = 
eMc^ and the expression relating Lq to the time-averaged 
mass accretion rate, Lq = 9 x 10^^ <M>_n ergs s"', one ob- 
tains (see also Wijnands et al. 2001) Lq = [foutburst/(foutburst + 
fq)](<Loutbui-st>/130). Here, fq is the average time the source 
spends in quiescence between outbursts, <Loutburst> is the 
average luminosity during an outburst, and foutbui-st is de- 
fined above. For XTE J2123-058, we used the RXTE/ASM 
light curve shown in Figure (SJl to determine that the mean 
1.5-12 keV count rate during the 40 days (foutburst) between 
MJD 50,990 and MJD 51,030 was 3.2 s"', corresponding to 
<Loutbui-st>= 1-1 X lO-'^ ergs s~' for li = 8.5 kpc. From the lu- 
minosity upper limit of Lq < 1.4 x lO-'^ ergs s"', we obtain 
a lower limit on the recurrence time, foutburst + fq, of 67 years. 
The known outburst history of XTE J2123-058, one outburst 
in 33 years, assuming that no outbursts were missed, is con- 
sistent with these long predicted recurrence times. However, 
if the recurrence time is in fact shorter than 67 years, then a 



mechanism of enhanced cooling of the core would be neces- 
sary. The most often mentioned mechanism is neutrino cool- 
ing of the core due to the direct Urea process, and this mech- 
anism r equires a ne utron star with a mass higher than 1.7- 
1.8M0(Colpi et al. 2001). The mass of the neutron star in 
XTE J2123-058 has been measured at 1.5 ±O.3M0 or 1.04- 
1.56M0, so that the mass would have to be at the upper end 
of the error range for the direct Urea process to operate. 

While we have been focusing primarily on the implica- 
tions of the upper limit on the temperature and luminosity 
of the thermal component, it is notable that the power-law 
fit to the spectrum of XTE J2123-058 also indicates that its 
X-ray luminosity, at (9;|;3) x 10-" ergs s"' in the unabsorbed 
0.3-8 keV band, is among the lowest for neutron star SXTs. 
The quiescent luminosity of SAX J1810.8-2609, « 1 x 10^^ 
ergs s"' ( Jonker. Wiinands & van der Klis 2004), is essen- 
tially the same as for XTE J2123-058, and the only other sys- 
tem in this class known to have a lower quiescent X-ray lumi- 
nosity, SAX J1808.4-3658, is also the only millisecond X-ray 
pulsar system for which there is a measurement of its quies- 
cent X-ray spectrum. The XMM-Newton spectrum of SAX 
J1808.4-3658 is well described by a power-law with F = 1.5 
and a 0.5-10 keV unabsorbed luminosity of 5 x 10^' ergs s~', 
and the upper limit on the contribution from a thermal com- 
ponent is 10% of the total X-ray luminosity (Camoan a et alJ 
2002). If the spectrum of XTE J2123-058 consists of only 
a power-law component, then its F = 3.l;|;og is quite a bit 
softer than that of SAX J1808.4-3658, and the difference 
would have to be explained. On the other hand, if the spec- 
trum for XTE J2123-058 consists of a thermal component and 
a power-law, then its power-law could have the same slope 
as SAX J1808.4-36 58. Concerning die S AX J1808.4-3658 
thermal component. iCampana et al.f i2002h conclude that en- 
hanced neutron star cooling is required to obtain a quiescent 
luminosity as low as the XMM-Newton upper limit implies. 
However, as mentioned above, such conclusions depend on 
the assumption that the recent outburst behavior is indicative 
of the average mass accretion rate over the past 10,000 years. 

In summary, our X-ray measurements of XTE J2123-058 
in quiescence show that this system is one of the faintest and 
least luminous of the neutron star SXTs. A comparison of the 
XTE J2123-058 X-ray properties in outburst and quiescence 
to that of four other neutron star SXTs that exhibit short out- 
bursts shows that the systems with a higher degree of outburst 
activity tend to have more luminous thermal components in 
quiescence, and this is in-line with predictions of the theory of 
deep crustal heating. For XTE J2123-058, the upper limit on 
the thermal luminosity is consistent with this theory without 
requiring enhanced neutron star cooling if the outburst recur- 
rence time is >67 years, which is consistent with the known 
outburst history of this source. 
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Table 1 

XTE J2123-058 X-RAY AND OPTICAL Observations 



Observatory 



UT Date 



Energy Band Exposure Time (s) 



BeppoSAX 2000 May 12-13 0.1-lOkeV 46,340° 
MDM 2000 July 24 Harris R 3 X 600 

Lick 2002 September 1 1 Kron-Cousins R 2 X 700 

Chandra 2002 November 13 0.3-8 keV 17,706 
MDM 2002 November 25 Harris R 8x600 

"This is the MECS exposure time. The LECS exposure time is 17,080 s. 



Table 2 
Spectral Fits 



Model 




r 




kThb 




Fit 




(10^" cm-2) 




(10 ergs cm ^ s 


') (eV) 


(10^2 ergs s 


') Quality'' 


Pl 
bb 

bb+pl 


7 n+40.0 
'■"-7.0 

0.0^0,0 
8.0^3- 


3 1+2** 
2 5+2' 


1 9+31.9 
'■^-0.5 

70+3'*' 

— 7"j.f 


260+-''0 


4Q+2.76 
"■^^-0.12 


0.69 
0.97 
0.66 


pl 
bb 

bb+pl 


6.0 
6.0 
6.0 


-'■'-0.6 
2 5+2-' 


'■'-0.4 
1^-^ 


240+11; 
' '"-90 


63+<'-22 
"■"-0.5 


0.69 
0.987 
0.64 



"Errors are 90% confidence for all parameters. 

''O.S-S keV unabsorbed flux. At d = 8.5 kpc, = lO"'"* ergs cm"^ g-l coiTesponds to 8.6 X 10^' ergs s"' . 
'^Bolometric luminosity at li = 8.5 kpc. 

''Fraction of 10,000 simulated spectra for which the C-statistic was better than the C-statistic obtained for the actual spectrum. 



Table 3 
NSA PLUS Power-law Fits 



d 

(kpc) 


(eV) (10" 


Fnsa" 

-14 - 
ergs cm s 


1 '' 
') (10^2 gj.gj, j,- 


r 


(10 ergs cm ~ s ') 


6 

8.5 
11 


<57 
<66 
< 73 


< 1.1 

< 1.1 

< 1.0 


< 0.80 

< 1.4 
<2.1 


7 ^+1.4 
^■-"-2.5 
2 5+'-'' 


56+0^2 
"•JO_0.42 

63+'"*2 

"■°-^-0.60 
77+0.65 
"■"-0.77 



"0.3-8 keV unabsorbed flux. 

''Bolometric luminosity as seen by a distant observer using 47rR2(j7j*^(l -2GM /Rc^), where R = 10 km and M = 1.4Mq. 
'^0.3-8 keV unabsorbed flux. Atd = 8.5 kpc, Fpi = lO"'"* ergs cm"- s"' corresponds to 8.6 X 10^' ergs s"'. 



Table 4 

NSA Parameters for Short Outburst Field NS SXTs and X-Ray Activity 



Source" 


d 




R 


Loo'' <M>'' 






(kpc) 


(eV) 


(km) 


(ergss"') (M© yr" 


') 


Aql X-1 


5 


94-108'' 


13.2 


(5.3-9.4) X 10" 0.828 1.3 x 10" 


10 


4U 1608-522 


3.6 


170 ±30 


9.4 


5.3 X 10'-^ 0.748 6.9 x 10" 


11 


Cen X-4 


1.2 


76±7 


12.9 


4.8 x 10'2 0.824 1.9 x 10" 


11 


SAX J 1810.8-2609 


4.9 


< 72 


10 


<2.0x 10" 0.765 1.7 X 10" 


12 


XTE J2 123-05 8 


8.5 


< 66 


10 


< 1.4x 10" 0.765 2.9 X 10" 


12 



"Parameters for sources other than XTE J2123-058 are from lRutledge et alJ lT999tl2()0lll2002all and LTonker. Wiinands & van der KUd|2()04l) . 

''Bolometric luminosity as seen by a distant observer calculated using 'i-KR? aT^gg^ . 
"g = ^X-lGMjRc^ assuming M = 1 .4 M© . 

''The time-averaged mass accretion rate over the past 33 years, estimated according to <M>= sL^-^N, where « = 1.1 X 10"-^ s^ cm"^ (see Equation 1). 
''The range of values found in three separate Chandra observations. 



